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ABSTRACT In contrast to the rapidly increasing knowledge on genome content and diversity of bacterial viruses, insights in intra- 
cellular phage development and its impact on bacterial physiology are very limited. We present a multifaceted study combining 
quantitative PCR (qPCR), microarray, RNA-seq, and two-dimensional gel electrophoresis (2D-GE), to obtain a global overview 
of alterations in DNA, RNA, and protein content in Pseudomonas aeruginosa PAOl cells upon infection with the strictly lytic 
phage LUZ19. Viral genome replication occurs in the second half of the phage infection cycle and coincides with degradation of 
the bacterial genome. At the RNA level, there is a sharp increase in viral mRNAs from 23 to 60% of all transcripts after 5 and 
15 min of infection, respectively. Although microarray analysis revealed a complex pattern of bacterial up- and downregulated 
genes, the accumulation of viral mRNA clearly coincides with a general breakdown of abundant bacterial transcripts. Two- 
dimensional gel electrophoretic analyses shows no bacterial protein degradation during phage infection, and seven stress-related 
bacterial proteins appear. Moreover, the two most abundantly expressed early and late-early phage proteins, LUZ19 gene prod- 
uct 13 (Gpl3) and Gp21, completely inhibit P. aeruginosa growth when expressed from a single-copy plasmid. Since Gpl3 en- 
codes a predicted GNAT acetyltransferase, this observation points at a crucial but yet unexplored level of posttranslational viral 
control during infection. 

IMPORTANCE Massive genome sequencing has led to important insights into the enormous genetic diversity of bacterial viruses 
(bacteriophages). However, for nearly all known phages, information on the impact of the phage infection on host physiology 
and intracellular phage development is scarce. This aspect of phage research should be revitalized, as phages evolved genes which 
can shut down or redirect bacterial processes in a very efficient way, which can be exploited towards antibacterial design. In this 
context, we initiated a study of the human opportunistic pathogen Pseudomonas aeruginosa under attack by one its most com- 
mon predators, the Phikmvlikevirus. By analyzing various stages of infection at different levels, this study uncovers new features 
of phage infection, representing a cornerstone for future studies on members of this phage genus. 
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The outcome of a virulent bacteriophage infection is deter- 
mined by a complex struggle between the virus and the bacte- 
rial host during all stages of infection (1). After adsorbing to a 
susceptible host, the phage ejects its DNA into the bacterium and 
attempts to establish a favorable climate for phage production. 
This process occurs worldwide approximately 10 23 times every 
second (2) and leads to a highly dynamic coevolutionary process 
of resistance and counterresistance in virus and host. Although 
this process is one of the major forces in microbial evolution (3), 
current knowledge of the impact of phage infection on host phys- 
iology is limited to a select number of (mainly) Escherichia coli 
phages. Well-known examples of these phage-host interactions 
are the transcriptional shutoff in E. coli by the Gp2 and Ale pro- 



teins of phages T7 and T4, respectively, or the cessation of DNA 
replication in Staphylococcus aureus through action of Gpl04 of 
phage 77 (4-6). 

This paper focuses on the intracellular development of the 
strictly lytic phage LUZ19 infecting Pseudomonas aeruginosa 
PAOl, an isolate of the Pseudomonas phage phiKMV (7). This 
podovirus carries a 43,548-bp genome, which is delineated by long 
direct terminal repeats. It contains 54 predicted genes organized 
in clusters of early, middle, and late expression. The genome is 
over 90% identical to other members of the Phikmvlikevirus ge- 
nus, which is a genus of extremely abundant Pseudomonas phages 
which can be readily isolated from wastewaters around the world 
(8). These phages belong to the family of the Autographivirinae, as 
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they encode their own single-subunit RNA polymerase (RNAP), 
relying on the transcriptional apparatus of the host only for tran- 
scription of the early and middle genes (9). 

The host organism in this study, P. aeruginosa PAOl, is a bac- 
terium with extreme metabolic versatility. It is able to colonize 
various environmental habitats and to persist and grow under 
nutrient-poor and hostile conditions (10). P. aeruginosa can 
switch from a soil and water organism to an opportunistic human 
pathogen and is responsible for life-threatening infections in cys- 
tic fibrosis and burn wound patients (11). The extraordinary 
physiological capacities of this Gram-negative bacterium have 
been studied intensively, mainly by using the available P. aerugi- 
nosa GeneChip (Affymetrix) to analyze differential gene expres- 
sion (e.g., see references 12 to 14). 

In this study, we combine quantitative PCR (qPCR), microar- 
ray, RNA-seq, and two-dimensional gel electrophoresis (2D-GE) 
proteome methods in the study of the global effect of LUZ19 in- 
fection on host physiology. By analyzing different stages of infec- 
tion, we gather insights into viral transcriptional overtake and 
stress responses provoked in the host and show the influence of a 
virulent phage infection on the DNA, RNA, and protein content of 
P. aeruginosa. 

RESULTS 

Characterization of LUZ19 infection. To analyze phage-host in- 
terplay at various stages of infection, we first studied the gross 
infection parameters of LUZ19. Exponentially growing P. aerugi- 
nosa cells were infected at an OD at 600 nm (OD 500 ) of 0.3 ( 1 .2 X 
10 8 CFU/ml). Stepwise increasing amounts of phages added to 
these cultures clearly showed that a single infection cycle is com- 
pleted in approximately 24 min and that addition of a 10-fold 
excess of phages leads to coordinated cell death after one genera- 
tion. Moreover, within the first minute after phage addition, a 
4-log reduction in nonadsorbed phage particles was observed (k a 
= 3.34 X 10~ 7 ml min -1 ), which is indicative of synchronized 
infected cultures (see Fig. SI in the supplemental material). For all 
further experiments, a 10-fold excess of phages over bacteria was 
used. 

Viral DNA amplification initiates after 14 min. Phage LUZ19 
encodes a replication cluster which closely resembles the well- 
studied phage T7 replisome and initiates genome amplification 
6min after the onset of infection (8, 15, 16). To study the timingof 
genome replication and the fate of P. aeruginosa DNA during 
LUZ19 infection, we extracted total DNA from infected cells at 5- 
to 7-min intervals. Next, we determined the abundance of 
genomic copies of gene 19 (phage) and three bacterial reporter 
genes, oprL, IppL, and PA3326, by quantitative PCR. In this exper- 
iment, we observed that LUZ19 genome replication starts only 
after 14 min of infection. Host DNA replication is not shut down 
before this time point, as the copy number of the selected genes 
doubles compared to that at the moment of infection (Fig. 1; see 
also Fig. S2 and S3 in the supplemental material). However, coin- 
ciding with phage genome amplification, the bacterial genome is 
degraded, presumably providing the phage with nucleotide build- 
ing blocks (Fig. 1A). As this degradation is observed for the three 
distantly carried bacterial reporter genes, it more than likely rep- 
resents a genome-wide phenomenon. Analogous to phage T7 
(17), this breakdown of the host chromosome is probably medi- 
ated by the combined action of the 5' to 3' exonuclease (gene 22) 
and the type 7 endonuclease (gene 23), which are encoded among 




Time after infection (min.) 

FIG 1 qPCR analysis of LUZ19-infected P. aeruginosa cells in LB medium (A) 
andM9 minimal medium (B). Cellular DNA was extracted at regular intervals, 
and the quantity of genes 19 (LUZ19, black line) and oprL (P. aeruginosa, gray 
line) over time was monitored in three independent reactions. The moment of 
lysis is indicated with an asterisk. The results obtained for two other bacterial 
reporter genes (UpL and PA3326) are shown in Fig. S3 in the supplemental 
material. 



the DNA replication genes of LUZ19. In LB medium, approxi- 
mately half of the P. aeruginosa PAO 1 genome is degraded at the 
end of the infection cycle. Interestingly, the entire host chromo- 
some is broken down at the moment of cell lysis when this exper- 
iment is performed in M9 minimal medium (Fig. IB). 

Transcript isolation and sequencing. In a next step, we ana- 
lyzed changes in the cellular RNA content provoked by phage 
infection. In two independent experiments, total cellular RNA was 
extracted immediately before and after 5 and 15 min of phage 
infection. rRNAs were depleted by subtractive hybridization, and 
mRNA samples were checked for integrity by an Agilent Bioana- 
lyzer. These samples were subjected to both classical microarray 
(Affymetrix) and RNA-seq analysis. In the latter analysis, between 
16 and 51 million reads of more than 20 bases were generated, 
resulting in the sequencing of 780 to 2,396 million bases in total 
(see Table S 1 in the supplemental material) . By iterative alignment 
and 3 '-end trimming, on average 88.5% of the reads could be 
mapped to either the bacterial or the phage genome. Although we 
performed subtractive hybridization to remove bacterial rRNAs as 
described in Materials and Methods (18), the majority of sequence 
reads (up to 98%) still represented rRNAs and other structural 
RNAs. This inability to successfully extract rRNAs from P. aerugi- 
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FIG 2 Differential expression of P. aeruginosa genes after LUZ19 infection. Depicted are the total number of genes which displayed a minimal 2-fold up- 
downregulation (green and red, respectively) in the microarray study, categorized according to the work of Stover et al. (21). 



nosa total RNA by this method was recently also reported by 
Dotschetal. (19). 

Influence of phage infection on bacterial gene expression. 

The microarray analysis showed that 220 P. aeruginosa genes are 
relatively upregulated during LUZ19 infection (posterior proba- 
bility of differential expression [PPDE], >0.95). This is distinctly 
different from other lytic phage-host systems where few, if any, 
host genes become activated (17, 20). Moreover, 13 operons were 
coordinately downregulated, also suggesting targeted shutdown 
of transcription (Fig. 2; see also Table S2 in the supplemental 
material). By classifying these genes in various functional groups 
(2 1), we observed that 25 significantly regulated genes are directly 
classified as energy metabolism genes. However, many more are 
indirectly linked to the energy status of the cell, such as enzymes 
and other macromolecule metabolism genes. Two of the most 
highly induced genes, glyAl (+4.97) and hutU (+4.02), for in- 
stance, play an important role in amino acid metabolism. Another 
striking observation is the rapid downregulation of the nuo gene 
cluster, encoding the NADH dehydrogenase complex, and the re- 
pression of the cytochrome o ubiquinol oxidase subunits. 

Also notable are the changes in gene expression of various 
membrane-related genes. For example, the wbp and am operons, 
which mediate the biosynthesis of the lipopolysaccharide (LPS) O 
antigen and lipid A layer, respectively, are strongly downregulated 
during phage infection. Moreover, three type IV pilus genes, 
whose products serve as primary receptors for LUZ19 (22), are 
repressed. We refer to the supplemental material for a detailed 
overview and discussion of all differentially expressed genes and 
operons. 

Instability of abundant bacterial transcripts. Two major ad- 
vantages of RNA-seq over the GeneChip array data are that it 
allows identification of phage transcripts and estimation of tran- 
script abundance. In comparing the total amounts of nonstruc- 
tural RNA reads originating from phage and bacteria from our 
RNA-seq data, we observed a steep increase from 23.6 to 60.2% of 



phage mRNA after 5 and 15 min of infection, respectively. As this 
increase could be caused by both rapid phage transcription and 
host mRNA degradation, we analyzed the RNA-seq data sets for 
differential gene expression with two algorithms, edgeR and 
DESeq (23, 24) at the default false discovery rate of 0.1 (25) 
(Fig. 3A). For the bulk of the bacterial transcripts, no obvious 
degradation was observed; 49 and 52% of them were at least as 
abundant after 5 and 15 min of infection, respectively, as in the 
uninfected sample. However, it is clear that the most abundant 
host transcripts are degraded when the phage infection proceeds, 
except for a few (conserved) hypothetical gene transcripts: 
PA4463, PA3623, PA4323, PA0257, and PA0689 (Fig. 3A). The 
abundant transcripts on the right end of the graph represent 
rRNA, which appears unaffected by LUZ19 infection. 

Phage transcription. Using a pile-up of all RNA sequence 
reads of phage transcripts, we visualized the genome coverage at 
two points during infection. This profile shows the expected shift 
in transcripts from early to late genes after 5 and 15 min of infec- 
tion, respectively, indicative of the temporal pattern of gene ex- 
pression (Fig. 3B). Moreover, a characteristic enrichment of se- 
quence reads in promoter regions is observed. This phenomenon 
arises as the 5' end of a transcript is always present in the leftmost 
fragment after random fragmentation. As such, independently of 
the exact fragmentation point, reads which originate from the 
transcript left ends are overrepresented (19). 

Early transcription is clearly initiated from three promoters at 
positions 913, 1147, and 1248 in front of the first gene. These 
promoters strongly resemble the bacterial cr 70 consensus sequence 
(TTGACA at the -35 box, TAcAAg at -10). The second viral 
transcriptional profile, generated after 15 min of infection, 
showed various sharp peaks throughout the genome. The major- 
ity of these correspond to intergenic region upstream genes 23, 27 
(downstream of the RNA polymerase gene), 32 (capsid protein), 
and 37, 42, and 47 (head decoration proteins). Although one 
could hypothesize that these transcripts originate from virus- 
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FIG 3 Analysis of the RNA-seq data. (A) Differential expression in P. aeruginosa genes by DESeq analysis of gene counts at a false discovery rate of 0. 1 . Each dot 
represents a specific gene after 5 (blue) and 15 (red) minutes of LUZ19 infection, with its estimated transcript abundance on the x axis and its log 2 fold change 
compared to the uninfected sample on they axis. The dots on the right end of the graph represent rRNAs. (B) Pile-up of LUZ19 reads after 5 (top graph) and 15 
(bottom graph) minutes of infection. The intergenic regions which are presumably associated with host- and phage-specific promoters are indicated with red and 
blue triangles, respectively. The locations of single-stranded interruptions on the noncoding strand as defined by Kulakov et al. (26) are shown as green triangles. 



specific RNA polymerase, no common promoter motif could be 
detected at these sites. The only detectable intergenic motif (CCT 
ACTCCGG) at positions 3132, 6755, 14415, and 24767 was previ- 
ously associated with localized, single-stranded nicks on the non- 
coding strand (26). However, no relation to transcript abundance 
is apparent (Fig. 3B). 

Subtle changes in bacterial proteome after phage LUZ19 in- 
fection. In a next step, we compared 2D-PAGE cytoplasmic pro- 
teome maps of exponentially growing, uninfected P. aeruginosa 
PAOl cells to corresponding maps of cells after infection with 
LUZ19. The reference map of uninfected P. aeruginosa cells con- 
tained (on average) 1,128 spots per 200 fxg of cell lysate and was 
recently published by our group (27). Using similar techniques, 
three phage infection maps were created at 5-min intervals, giving 
a spread of analysis points during the entire phage cycle (see 
Fig. S4 in the supplemental material). By comparing the 2D maps, 
it is apparent that the global impact of phage infection on the host 
proteome is limited. With only one exception (see below), we did 
not detect protein modifications that alter the protein mass or 
isoelectric point. Moreover, almost no degradation of host pro- 
teins was observed. In the pi range of 4 to 7, 12 new spots were 
detected during infection and identified using electrospray 
ionization-tandem mass spectrometry (ESI-MS-MS). Six of these 
represented host proteins are involved in stress responses (PasP 
protease, PA2807, and AotJ), energy metabolism (SdhB), and 



translation (EF-Tu and RplL). An intriguing observation con- 
cerns RpoA (PA4238), the alpha subunit of the bacterial RNA 
polymerase (RNAP). In the noninfected sample, two isoforms of 
this protein were detected, both with an estimated protein mass of 
44 kDa. LUZ 1 9 infection induced a shift of this spot to a molecular 
mass of approximately 32 kDa, indicative of degradation upon 
phage infection or sample preparation. We refer to the supple- 
mental material for a detailed discussion of these proteins. 

Appearance of abundant phage proteins. At distinct time 
points during infection, five LUZ 19 proteins were detected on the 
two-dimensional electrophoresis (2-DE) infection maps. They 
belong to different gene classes (early, middle, and late) and again 
illustrate the timed regulation of phage infection. First, a highly 
abundant early protein (GplO) appeared (Fig. 4, inset). Shortly 
thereafter, two late-early/middle proteins, Gpl3 and Gp21, were 
detected. Based on the location of the protein spot in the 2D gel 
and identified peptides during mass spectrometry, the gene anno- 
tation of Gpl3 was reevaluated (see the supplemental material). 
The most abundant phage protein was, as expected, the major 
capsid protein Gp32, which appeared first after 10 min of infec- 
tion and increased in abundance in the course of the phage infec- 
tion. The putative tail fiber Gp41 was also identified. Other phage 
proteins either were too small to be detected, were not abundant 
enough, or fell outside our pi range. 

Based on sequence similarity, the function of GplO and Gp21 
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FIG 4 Two-dimensional proteome map of P. aeruginosa cells, 5 min after infection with LUZ19. The right corner shows outtakes of two corresponding gel 
regions from the uninfected sample (upper panel) and after 5 and 15 min of infection, respectively (lower panel). The encircled spots were identified as LUZ19 
GplO (left) and the major phage capsid protein Gp32 (right). 



cannot be predicted, although the latter is conserved in many 
podoviruses and prophages (PHA02030, E = 2.1 X e -151 ). In 
contrast, Gpl3 is highly similar to various bacterial GNAT acetyl- 
transferases (pfaml3508, E = 2.51 X e~ 8 ). To gain insight into the 
individual impact of these abundant early phage proteins on host 
physiology, they were cloned in a single-copy pUC18-mini-Tn7T- 
LAC expression vector (28) which integrates into the bacterial 
genome (Fig. 5A), and cell growth was monitored upon induc- 
tion. Interestingly, while GplO had no influence on P. aeruginosa 
growth, expression of both Gpl3 and Gp21 was sufficient to block 
bacterial growth completely in both solid and liquid media 
(Fig. 5B). Time-lapse microscopic recordings of these cells 
showed that the presence of the phage acetyltransferase slows bac- 
terial growth down, while Gp21 expression leads to an immediate 
growth arrest (Fig. 5C; see also movies Ml to M3 in the supple- 
mental material). Given that Gp21 is encoded among genes in- 
volved in DNA replication and directly precedes the phage- 
encoded exo- and endonucleases, we explored the possibility that 
this protein is involved in shutdown of bacterial DNA replication. 
Therefore, we followed the incorporation of radiolabeled macro- 
molecular precursors in P. aeruginosa cells with and without in- 
duction. For Gp21, we observed an immediate drop of 20 to 25% 
in transcription and translation. However, DNA replication 
seemed unaffected (Fig. 5D). In contrast, translation appears un- 
affected in the first 30 min after Gpl3 induction, while DNA rep- 
lication and transcription levels decrease approximately 30% 
compared to those in noninduced P. aeruginosa cells. Only 21% 
residual transcriptional activity is remaining after 65 min of Gpl3 
induction (Fig. 5D). The elucidation of the binding partner(s) and 
working mechanisms of these proteins is under investigation. 



DISCUSSION 

In this paper, we describe a study of the total DNA, RNA, and 
cytoplasmic protein content of phage-infected P. aeruginosa 
PAOl cells based on qPCR analyses, whole-cell proteomics, and 
transcriptome analyses by both microarray and RNA-seq technol- 
ogy. We selected the obligatory lytic podovirus LUZ19 for this 
analysis, which adsorbs rapidly to the host cell and completes its 
infection cycle within 25 min. Early bacteriophage transcription 
could be clearly associated with three promoters directly upstream 
of gene 1, which closely resembles the a 70 consensus sequence. 
Later in infection, many potential transcriptional start sites could 
be detected alongside the phage genome, but no LUZ19-specific 
promoter consensus sequence could be retrieved. Due to the 
method of cDNA library construction (see the supplemental ma- 
terial), no information on antisense transcripts could be inferred. 
In the future, directional transcriptome maps for this (and other) 
phages can be generated by constructing libraries from first- 
strand cDNA(29). 

A clear observation was the rapid accumulation of phage tran- 
scripts during LUZ19 infection, in contrast to breakdown of bac- 
terial mRNA. Experiments presented in this study clearly show 
that this breakdown initiates before bacterial genome replication 
arrests. There is certainly an advantage for a virulent phage to 
accelerate host mRNA degradation immediately after infection, as 
this contributes to making the translation apparatus more rapidly 
available for viral mRNAs and facilitates the transition from the 
host to phage gene expression (30). In the related coliphage T7, the 
protein kinase Gp0.7 phosphorylates the C-terminal domain of 
RNase E, which protects T7 mRNA from degradation (31). As no 
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FIG 5 Analysis of phenotypic effects provoked by LUZ19 genes 10, 13, and 21 on P. aeruginosa PAOl. (A) The three genes were stably incorporated in single 
copy in the bacterial genome through Tn7-mediated integration of the expression vector pUC18-mini-Tn7T-LAC. ORF, open reading frame. (B) Phenotypic 
effect of phage gene expression on bacterial growth on solid (left) and in liquid (right) (gray line, empty vector; squares, Gpl3; dots, Gp21) LB media 
supplemented with 1 mM isopropyl-j3-D-thiogalactopyranoside (IPTG). (C) Time-lapse screen shots of a bacterial culture at the moment of gene induction 
(upper pictures) and after 6 h of induction at 37°C (bottom pictures). (D) Incorporation of radioactively labeled [ 3 H] uridine, [ 3 H]leucine, and [ 3 H] thymidine, 
followed in time upon induction of Gpl3 (left panel) and Gp21 (right panel). All experiments were performed in triplicate; standard deviations are indicated on 
the graphs. 



homologue to the T7 kinase is encoded in LUZ19, we speculate on 
another type of modulation of the host degradosome. 

It is possible that LUZ19 exerts posttranslational control by 
acetylation instead of phosphorylation. The early phage protein 
Gpl3 is an acetyltransferase for which we showed that single-copy 



expression is detrimental for P. aeruginosa and hinders both DNA 
replication and transcription. These types of enzymes are recog- 
nized more and more as serving pivotal roles in the posttransla- 
tional control of central microbial metabolism and are reported to 
be able to modulate enzyme activity (32, 33). Phage-induced acet- 
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ylation of specific phage or host proteins has not been reported 
thus far, and the study of the targets of Gp 13 is a very interesting 
area of future research. 

Even though virulent phage infection has a clear impact on the 
cellular RNA content, only moderate changes were observed at the 
protein level. This is not unexpected, as the average half-life of 
bacterial proteins exceeds the time span of the phage infection 
(34). Possible host transcription shutdown or bacterial mRNA 
degradation would in this case not result in a decrease of protein 
spots on 2-DE gels. We observed significant upregulation of a 
small number of proteins which previously have been shown to be 
associated with stress, like the EF-Tu elongation factor and the 50S 
ribosomal protein RplL (35, 36). In addition, a notable result from 
this study was the observed mass shift in the alpha subunit of the 
RNA polymerase during phage infection. All bacteriophages 
known so far redirect the bacterial RNAP toward transcription of 
the viral genome. While the majority of them recruit the host 
RNAP throughout the infection cycle, some phages (like LUZ19, 
T7, and N4) encode their own RNAP and rely on the host's ma- 
chinery only for transcription of early or late genes (17, 37). For 
this molecular hijacking, two main types of phage-induced mod- 
ifications have been described so far: covalent RNAP modifica- 
tions and modifications through RNAP -binding proteins (38). 
This observation could represent an entirely novel mechanism of 
viral control over the host transcriptional apparatus. 

MATERIALS AND METHODS 

Unless stated otherwise, Pseudomonas aeruginosa PAOl cells grown in LB 
medium were used for all experiments. Inducible single-copy expression 
of phage genes was achieved by integration of the target gene under a lac 
promoter into the bacterial genome using the pUC18-mini-Tn7T-LAC 
and pTNS2 suicide plasmids (28). High-titer phage stocks of LUZ19 were 
purified using two consecutive CsCl gradient centrifugations. 

For the transcriptome analyses, total RNA from phage-infected cells 
was isolated from the bacteria by a modified hot-phenol method in com- 
bination with the RNeasy midikit (Qiagen, Hilden, Germany). rRNA was 
extracted by hybridization to antisense oligonucleotides, followed by pull- 
down through binding by magnetic beads (MicrobExpress kit; Ambion). 
Subsequent cDNA library constructions for microarray and RNA-seq 
analysis were based on standard protocols for the Affymetrix Pseudomo- 
nas GeneChip and Illumina HISeq2000 sequencing, respectively. Array 
analyses were carried out using R and Bioconductor, and the preprocessed 
expression data were analyzed by a regularized t test based on a Bayesian 
statistical framework using the CyberT algorithm. The raw sequence data 
were checked for quality and mapped to the P. aeruginosa PAOl 
(NC_002516.2) and LUZ19 (NC_010326.1) reference genomes using 
Bowtie (39). Count data were analyzed for differential gene expression 
using the DEB web interface. 

The proteome analysis was performed by halting the LUZ19 infection 
at 5-min intervals and performing two-dimensional electrophoresis (2- 
DE) separations and image analyses on clear cell lysates. 2-DE gels were 
silver stained, and image acquisition was performed using a calibrated 
flatbed ImageScanner, combined with LabScan software according to the 
manufacturer's protocols. 2-DE maps were analyzed, and spot data were 
generated by standard spot detection parameters using ImageMaster 2D 
platinum software. Differential spots were picked and identified using 
ESI-MS-MS. 

Detailed procedures can be found in the supplemental material. 
SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi: 10.11 28/mBio.0006 1 - 1 3/-/DCSupplemental. 
Text SI, DOCX file, 0.1 MB. 



Figure SI, PDF file, 0.2 MB. 
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